(N 



o 



in 



Optical Spectral Singularities and Coherent Perfect Absorption in a 

Two- Layer Spherical Medium 

Ali Mostafazadeh* and Mustafa Sansamant 
Department of Mathematics, Koc University, 



^ | 34450 Sariyer, Istanbul, Turkey 

Abstract 

. An optical spectral singularity is a zero-width resonance that corresponds to lasing at threshold 

gain. Its time-reversal causes coherent perfect absorption of light and forms the theoretical basis of 
antilasing. In this article we explore optical spectral singularities of a two-layer spherical medium. 
. In particular, we examine the cases that a gain medium is coated by a thin layer of high-refractive 



index glass and a spherical glass covered by a layer of gain material. In the former case, the coating 
reduces the minimum radius required for exciting spectral singularities and gives rise to the formation 
of clusters of spectral singularities separated by wide spectral gaps. In the latter case, the coating 



CZ2 

Q " leads to a doubling of the number of spectral singularities. 

• i— i ' 

PACS numbers: 03.65.Nk, 42.25.Bs, 42.60.Da, 24.30. Gd 

'■^ ■ Keywords: Complex potential, spectral singularity, zero- width resonance, gain medium, coated spher- 

ical dye laser, coherent perfect absorption, antilaser 

^ ■ 1 Introduction 

The discovery pQ that the mathematical concept of a spectral singularity [2] has physical realizations 
as zero-width resonances of complex scattering potentials has motivated a detailed study of this phe- 



■ nomenon [3] - [12]. In particular, it is shown that optical spectral singularities (OSS) correspond to the 
CN ■ lasing at the threshold gain [8] and that a time-reversed optical spectral singularity [7] yields a coherent 
T 7"! . perfect absorption (CPA) of light, i.e., an antilaser [13] . 

. . Typical lasers are photonic devices consisting of an active medium placed inside an optical cavity. 

A particularly interesting type of lasers are those based on spherical granules where the surface of the 
sphere acts as the cavity [T3] . These are characterized by their extremely high quality factors and small 



volumes of excitation. Owing to these properties, a spherical laser can also be realized in which the 
active medium is located outside an spherical core [15]. 

In Ref. [11] we studied the OSS of a uniform spherical gain medium and showed for the radial 
(transverse) modes of a concrete spherical dye laser that the emergence of an OSS puts a lower bound 
on the radius of the gain medium. This is the minimum radius a*- mm ) required for lasing in these modes. 
In the present article we examine OSS and CPA for an active medium consisting of a spherical inner 
core and a spherical outer shell with different refractive indices, as shown in Figure Q3 In particular, 
we wish to explore the prospects of reducing the value of a*- mm ) by means of coating the spherical gain 
medium by a material with higher refractive indexQ 



*E-mail address: amostafazadeh@ku.edu.tr, Phone: +90 212 338 1462, Fax: +90 212 338 1559 

tE-mail address: msarisaman@ku.edu.tr, Phone: +90 212 338 1378, Fax: +90 212 338 1559 

1 It is well-known that the presence of such a coating increases the system's quality factor [151 116| . 
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Figure 1: (Color online) Two-layer spherical gain medium. a\ and 02 are the inner and outer radii of 
the spherical shell; n 1 and ri2 are the complex refractive indices of the inner core and the author shell, 
respectively. 



2 Radial Transverse Spherical Electromagnetic Waves 

Consider an optically active material with an inner spherical core of radius a\ and an outer spherical 
shell of thickness 02 — a\ placed in vacuum. Let n 1 and xi2 denote the complex refractive indices of the 
inner core and the outer shell, respectively, and suppose that they are independent of space and time. 
The electromagnetic (EM) waves interacting with this system satisfy the Maxwell equations: 

V-l5 = 0, V-B = 0, (1) 

VxE + B = 0, V x H - D = 0, (2) 

where D := £oi(r)E, H := ^^B, £q and (jlq are respectively the permeability and permittivity of the 
vacuum, r := \r\ is the radial spherical coordinate, 



l{r) :- 



nj for r < a\, 

tt| for «i < r < 02, (3) 
1 for r > 0,2, 



and each over-dot represents a time-derivative. According to ([2]), the electric field E = E(r,t) is a 
solution of the wave equation: 



E(r,t) + SfE(r,t) = 0, 01,02, (4) 

where Vt 2 := c 2 3(r) _1 V x Vx and c = (eo/^o)" 1 ^ 2 is the speed of light in vacuum. 

For a time-harmonic EM field with angular frequency uj that propagates in a charge-free medium, 
we have E(f,t) = e~ tLOt E(f), and (j4|) reduces to the time-independent Schrodinger equation: 

- X7 2 E(r) + v{r)E(r) = k 2 E(r), (5) 

where k := uj/c is the wave number and v is the complex barrier potential: v(r) := k 2 [l — 3(r)]H 

Following the analysis of Ref. [11] , we investigate transverse radially propagating spherical solutions 
of ([5]) that have the form 

E(f) = E(r)4>. (6) 



2 Note that ((5]) is nothing but the well-known Helmholtz equation, and the potential appearing in it is energy-dependent. 
For the purposes of our investigation this does not cause any difficulty. 
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Here <j> is the unit vector associated with the azimuthal angular coordinate eft of the spherical coordinate 
system. Inserting ([6]) in ([5]) yields 



— + -— + k 2 

dr 2 r dr 



v(r) - — 



E(r) = 0. 



(7) 



For r < a\, a\ < r < a 2 and r > a 2 where v takes constant values, we can transform ([7]) to the spherical 
Bessel equation of order u := [T7]. Therefore, 



E(r) 



A\j v (k\r) + B\n v (k\r) for r < a\, 
A 2 j v {k 2 r) + B 2 n u (k 2 r) for ai < r < a 2 , 
A 3 h { J-\kr)+B 3 h^(kr) for r > a 2 , 

,0 



(8) 



where the and 5j are constant numerical coefficients, j u ,n u and hi, are respectively the spherical 
Bessel, Neumann, and Hankel functions, and fcj := rij k for i = 1,2. Notice that the condition that the 
electric field be regular at the origin implies that B\ = 0. Consequently, E(0) = 0. 

Having obtained the explicit form of the electric field, we can compute the magnetic field using ([2]). 
This gives 



B(f,t) 



-ILOt 



B(r)8, 



where 6 is the unit vector associated with the spherical polar coordinate 



B(r) := iuj~ x E(r) = < 



i UJ 



l UJ 



[Axj v {kxr) + Bxn^kxr)] for r < ai, 
[A 2 j u (k2r) + B 2 h v (k 2 r)\ for a\ < r < a 2 , 



A 3 h { 1 } ) '(At) + B 3 h¥> \kr 



(2), 



(9) 



(10) 



for 



r > a 2 , 



and for each differentiable function / we define / according to /(r) := + ^) f(r). Notice than for 
every pair of differentiable functions /, g : M + —> R and positive real numbers r, s G R + , 



/(r)ff( S ) - g(s)f(r) = f(r)g'(s) - g(s)f'(r). 



(11) 



Because this quantity will frequently appear in our calculations, we denote it by f(r) d g(s) for brevity. 
In other words, 

(12) 
(13) 



f(r) d g(s) := f{r)g'(s) - g{s)f'{r) = f(r)g(s) - g(s)f(r) 



Recall that the Wronskian of / and g is given by W[f(r),g(r)] := f(r) d g{r). Therefore 

f{r)g{r) - g{r)f(r) = W[f(r),g(r)}. 



In order to relate the coefficients Aj and Bj appearing in (jHJ), we need to impose the appropriate 
matching conditions at the boundaries r = a\ and r = a 2 , [T7]- For the system we consider, these 
correspond to the condition that the parallel component of both the electric and magnetic fields must 
be continuous at the boundaries. In view of ([6]) and 0, this means that E = E(r) and B = B{r) must 
be continuous functions. We can satisfy this condition provided that we select the coefficients Ai and 
Bi such that E and B are continuous at r = a\ and r = a 2 . This gives 



Kn 

where for all p, q = 1, 2, 
Kp 9 := 



" A 1 ' 


= Ki2 


" A 2 ' 


K 22 


" A 2 ' 




" ^3 " 






= L 







B 2 


B 2 




B 3 



jvikqdp) n u{kqCLp) 
jv(kqQ"p) fii/(kqa,p) 



hP(ka 2 ) h { v\ka 2 ) 
~h£\ka 2 ) ~h { v\ka 2 ) 



(14) 



(15) 
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According to these equations and (|13p . 



det(Kpg) 
det(L) 



W\j„{x),n v (x)] 
W[h£Xx),h®{x)\ 



v q u p ) 



x=ka,2 



-i{ka 2 )~ 2 + 0. 



(16) 
(17) 



The fact that these quantities do not vanish was to be expected, because {jv,n v ) and (hP~\ h^) are pairs 
of linearly-independent solutions of a second order homogeneous linear differential equation. Eqs. (|16p 
and (I17p imply that K p(? and L are invertible matrices. We can use their inverse together with Eq. ([1 



and the fact that Ai ^ 0, to express A% and B3 in terms of A\ according to ^3 
M 2 \A\, where My are the entries of the transfer matrix 



MnAi and B 3 



M := L~ 1 K 22 K 12 1 Kii. 



(18) 



3 OSS and CPA for Radial Transverse Spherical Waves 



As discussed in Ref . [TT] , we can easily exploit the asymptotic properties of the spherical Hankel functions 
to infer that the reflection amplitude of our system is given by 



R :-- 



M 



21 



(19) 



Therefore, in order for characterizing OSS and CPA that correspond to the real poles and zeros of R, 
we only need to calculate Mu and M21. Using (fTTj) . ([TBI) . (fTTl) . (fT2j) . (fT5|) . (fT8j) . and doing the necessary 
algebra, we obtain 

M n = N 2 , M 2 i = -Nx, (20) 



where for both £ = 1,2, 
N e :=P 



n w (k 2 a 2 ) d h$ (ka 2 ) +Q jw(k 2 a 2 ) d h$ (ka 2 ) 



(21) 



P := i(kk 2 aia 2 ) j u {k 2 ai) d j v {k\ai), Q := i(kk 2 aia 2 ) j u {hai) d n v {k 2 a\). (22) 

^ y ^ y 

In view of (I20p - (I22p and the fact that f d g = —g d f, we can express (|19j) in the form 



[j'i/(fc2ai)^Jf (^lOi)] [^(fe^V^^c^)] - [/^(foai) Vj^fciOi)] [i„(fc2a 2 ) V^ J (fca 2 )] 



^(2), 



[ji/(fe20i) > 5'j 1/ (A; 1 oi)] [n„(/c 2 a 2 ) 3 K (ka 2 )] - [n^kia^d j^kiai)] [^(^2^2) V/ii, lj (fca 2 )] 



B 



Eqs. (|20p provide a simple demonstration of the fact that CPA corresponds to an OSS of the 
time-reversed system. To see this we recall that we can obtain the time-reversed system by complex- 
conjugating the refractive indices ri£. This implies 



ki — » k 



I-* 
Mn 



P -> P* 



-M 



21- 



<2 -> Q*, 
M 2 i -> 



-M 



in 



(23) 
(24) 



where we have used <^B), fl22j), (jUJ), (JTaJ), and the fact that /^(fco^)* = tiv'(ka 2 ). According to the 
last relation in (|24p . a CPA, that corresponds to R = 0, appears if and only if the reflection coefficient 
of the time-reversed system diverges, i.e., the latter develops an OSS and begins lasing at the threshold 
gain . Therefore, in the following we only consider the problem of locating OSS. We can easily obtain 



,(2), 
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the values of the physical parameters leading to a CPA of the spherical waves we consider by complex- 
conjugating the refractive indices ni and n 2 or changing the sign of the gain/attenuation coefficients 
of both the layers. 

A straightforward consequence of this observation is the fact that because the law of energy conser- 
vation prohibits the emergence of an OSS for the case that both the interior core and the outer shell of 
our system consist of lossy material, a CPA cannot be realized unless either the core or the outer shell 
includes a lossy medium. Furthermore, it is possible to generate both a CPA and an OSS, if our system 
involves both lossy and gain media. 

In order to determine the location of OSS in the space of the physical parameters of the system, we 
study the real zeros of M 2 i or alternatively N\ in the complex /c-plane. In view of (12ip and (122|) . this is 
equivalent to finding the real values of k fulfilling 

[jvfaaiftf j v (kiax)] [n l/ (k 2 ai)^ (ka 2 )] = [n u {k 2 ai)^ j u {kia 1 )][j u {k 2 ai)^ h^\ka 2 )]. (25) 

Noting that ki = n ik, this is a complex transcendental equation involving two complex variables, namely 
ni and n 2 , and two (positive) real variables: x\ := ka\ and x 2 := ka 2 . 

For a\ = a 2 and n 1 = n 2 that corresponds to a homogeneous spherical medium that we consider in 
|llj , the first factor on the left-hand side of (|25f) vanished identically and the first factor on its right-hand 
side becomes W[n u (x),j u (x)] with x = k±a± = x\n\. Because the latter is nonzero, (f25j) reduces to 

j u (k l a l fih^(ka l ) = Q. (26) 

We can use the identity: 

= ^ , (27) 



to express the derivative of the spherical Bessel, Neumann, and Hankel functions. Doing this in (|26|) 
gives rise to the equation for the spectral singularities of a spherical gain medium that we derive in |llj . 
namely 

d_ 

=<u dr 

Employing the identity (|27p in (|25p gives a more lengthy equation for OSS that we use in our 
numerical and graphical investigations. Before reporting the results of this investigation, however, we 
will carry out a perturbative analysis of (|25p . Similarly to the single-layer spherical medium we studied 
in [11], this turns out to reveal some basic properties of the solutions. 



^L\nhW(kr) = ^-lnj u (k ir ) . (28) 

ar r=ai dr r=a\ 



4 Perturbative Analysis of OSS 



Consider the Mie regime where a 2 > a\ 3> 2ir/k =: A and |rij| < 4. Then x 2 > x\ ^> 1, \kia\\ = 
|rii|xi = 27r|rii|ai/A 3> 1, and we can perform a large-xi (and -x 2 ) expansion of the terms appearing 
in (|25p . This requires using the following asymptotic expansions of the spherical Bessel, Neumann, and 
Hankel functions. 



ra„(z) 
^)(z) 



sm(z - ^ 



E 



l) s A 2s (u) cos(z-f ) ~ {-l) s A 2s+l (u) 

s=0 



?2s+l 



cos(z - f ) ~ {-l) s A 2s {v) , sin(z - f ) ^ {-l) s A 2s+ i{v) 



/2s 



+ 



?2s+l 



- f (_l)<(Z-Sl) 



^(-) s A 2s {v) 



-!) £ E 



/2s 



E 

s=0 



(-l) s A 2s+l {u) 



7 2s+l 



s=0 



s=0 



(29) 
(30) 
(31) 
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where 



r(v + k + i) _UtoH^ + k 



A k (y) :-- 



2 k k\ 



2 k k\T(v-k + l) ' 

and r stands for the Gamma function. 

Substituting (|29p - f)31 1) in (|25p . neglecting the quadratic and higher order terms in x^ 1 and x^ 1 in 
the resulting equation, noting that Aq(v) = 1 and A\(v) 



^Y^-, and introducing 



il := tan(/c^ a\ 



F(xi, m, n 2 ) := 



— ) = tan(xm£ - — ), 

m(n 2 -Ht 2 ) + n 2 ti(n 2 t 2 - i) 
n i (« — n 2 t 2 ) + n 2 ti(n 2 + zt 2 ) ' 



we find 



tan(x 2 n; 



~2' 



ni, n 2 ), 



(32) 
(33) 

(34) 



where we use to indicate that this equation is obtained by employing first-order perturbation theory. 
Solving ([Ml) for x 2 yields 



vri/ + 2tan- 1 (J r ) 1 



2n: 



2n. 



7r (2m + v) — i In 



1 + iJ r 
1-iF 



(35) 



where we have used "In" to denote the principal part of natural logarithm of its argument, suppressed 
the argument of T for brevity, and employed the identity 



tan 1 (z)=irm-\ In 

v ; 2% 



1 + iz 



1 



IZ 



m 



0,±1,±2, 



(36) 



The parameter m appearing in (|35[) is an integer that we identify with a mode number labeling OSS. 

For a homogeneous spherical medium, where a\ = a 2 and n i = n 2 , we have x\ = x 2 , T = — in i, 
and (|35|) gives 



2n 



7r (2m + v + l) — i In 



n x + 1 

m-i 



(37) 



This is in complete agreement with Eq. (23) of Ref. 

Because the left-hand side of (|35|) is real, we can express this complex equation as a pair of real 
equations: 



X2 — Re 
Im 

n 2 



1 



2n 2 

7r(2m + v 



ix(2m + u) — i In 



l + iJF 
1-iF 



0. 



i In 



1 + 



0. 



We also note that 



1 - iT , 

1 + (n 2 ti + im)(l + it 2 )(n 2 + l) 



(38) 
(39) 

(40) 



(n 2 ti - ini)(l - it 2 )(n 2 - 1) ' 

Next, we consider the special case of a coated spherical active medium with a 2 — a\ <S a±. In this 
case the gain/absorption properties of the coating can be neglected and n 2 may be assumed to take a 
real value that we label by ?t, 2 . This implies the equivalence of (I39j) and the condition that the absolute 
value of (|40p must be unity. Imposing this condition and noting that in this case t 2 is also real, we find 



n 2 ti + i n i 
n 2 ti — i n i 



n 2 



II 



n 2 + 1 



(41) 
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This is a real equation involving a complex variable, n i, and two real variables, n% and x\. In particular, 
it does not involve the mode number m. 

In order to use (|4ip for locating OSS, we express rti and ti in terms of their real and imaginary 
parts. Let rj\ and K\ denote the real and imaginary parts of n i , so that 

n 1 =r] 1 +inx, (42) 

and introduce 

a := xxr]x - — = kaxrjx - — = 2ir -J , (J := xx_K\ = kaxKx = — - — . (43) 

Then substituting ([4*2]) in ([32]) gives 

. . „. tana + itanh/3 . 

tx = tan a + tfi = — . (44) 

\ — i tan a tann p 

Next, we use (I42p . (|44p . and various trigonometric and hyperbolic identities to obtain the following 
explicit form of ([411) . 



/ 2 2 2 \ X 2 2 2\ 

r/i sinh(2/3) - ki sin(2a) + ( ??1 + + " 2 J cosh(2/3) + ( ^ + ~ " 2 J cos(2q) « 0. (45) 

V "2 + 1 / V "2 + 1 / 

A similar analysis reveals the fact that the following quantity is an integer. 

_! /l + iF\ _ 

m : = 7 - vr arg I - — — I , (46) 

where 

7 := — ^— + £ - (47) 

> -l , / 2n 2 [? ?1 sin(2a) + Kl sinh(2/3)] 1 

5 . 7T arCTan | [n 2_ ( ^2 + K 2 )]cosh(2 ^_ [n 2 + ^ 2 + K 2 ]cos ^ ^ 

and "arg(z)" denotes the principal argument of z that takes values in (— ir,ir]. The latter implies that 
m is one of the two integers satisfying the condition: 

7-l<77i<7 + l. (49) 

A more important implication of (|46|) is that it leads to the following explicit form of 



X(2m — m) 

a 2 -ai^a -\ , (50) 

4rt2 

where 

ao := £. (51) 

Note that the "arctan" appearing in (|48p stands for the principal value of "tan -1 " that takes values in 
[— §, §]• This in particular implies that |£| < 1/2. Therefore, 

Kl < (52) 
8n 2 

Another useful relation that follows from ([4"7|). ([58]) . and ([5T]) is 

4n 2 (ai+a ) . . 

7 = t v. (53) 
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Eq. (|50p is quite remarkable, for it indicates that if we choose ?7i,Ki,7i2 and x\ = 2irai/\ so that 
ii > 1 and ([4"5]) holds, then an OSS arises for a discrete set of values of the thickness of the coating. 
Because 2m — rh is an integer, according to ([50]) and ([521 . 



2m <J m > 7 — 1, «2 — a i ~ • (54) 

~ ~ 8ri2 

This means that the mode number m and the thickness a2 ~~ a i are bounded from below by (7 — l)/2 
and A/(8ri2), respectively. Notice that these bounds only depend on ai, A and 122- 

Next, we examine the physical implications of (|45|) for a typical optically active material that satisfies 



\ki\ < |/3| <C 771 < a. (55) 

In this case, we can ignore the terms of order two and higher in k\ and (3 in our calculations. Imple- 
menting this approximation in (|45p . using (|43p . and recalling that the imaginary part of ni, i.e., k±, is 
related to the gain coefficient g via k± = —gX/(4Tr), we find 



a l9 * 2 Z 7\~ 7^ U ° V ^ . (56) 



rjl + n\ + (rjf - ra 2 ,) cos(2a) 

mini + 1) 

This relation shows that the radius a\ of the inner core is inversely proportional to the gain coefficient. 
Furthermore, it puts curious upper and lower bounds on the possible values of a\ that are only sensitive 
to the gain coefficient and the real part of the refractive indices of the inner core and the outer shell: 

2r? 2 In 2 

"'mm < _ < max /^\ 

— ^ a\ ^ —5 — , {01) 

mg( n 2 + !) ^15(^2 + !) 

where n max and n m ; n are respectively the largest and smallest of r]i and 712. 

In view of (|57|) . the larger ri2 is, the smaller the lower bound of a\ gets. This confirms our expectation 
that coating a spherical gain medium by a material with higher refractive index reduces the lower bound 
on the radius of the gain medium. For example, if we take rti = 2.5 and choose the inner core to be 
made of a dye gain material with rji ~ 1.48 and g ~ 5 cm -1 , we find 

0.816 mm < a x S 2.330 mm. (58) 



Finally, we explore the consequences of (|55p . Implementing the above approximation scheme of 
neglecting second and higher order terms in f3 and k\ in (|50p yields 



«0 « TZZ- arctan \ Z2—T2 7Z2 , Z2WI77T3 ? ■ ( 59 ) 



A f 2r/ifi2 sin(2a) 

47rre2 \ rig — rjl — (r) 2 + n|) cos(2a) 

Taking a% = 1 mm to comply with (|58p . choosing A = 549 nm, r}\ ~ 1.48, and n2 = 2.5 as above, and 
using 0, g?D, (09]), H5K and (JH2J) we find 

a 2 - ai w [55 (2m - rh) + 27] nm, m w 7 = 18122. (60) 

The smallest allowed value of the thickness is therefore 02 — a\ ~ ao = 27 nm. 

In practice the thickness of the coating has a fixed value, and (j59[) and (|60p determine an approximate 
value of the mode number m: 

2n 2 (a 2 - ai - a ) , m , . 

mw Y-. (61) 

If we further approximate m by 7 and use the expression (|53p for the latter, then (|6ip gives 

m«^--. (62) 
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For example, for a 2 — a\ = 5 /an, this relation gives m ~ 9152, which coincides with the numerical result 
obtained directly from (|38j) . 

We can also use (|62f) to express the wavelength of the OSS in terms of the mode number. This gives 

A « p»L. (63) 
2m + v K ' 

The analogous expression for the case that we remove the coating is given by Eq. (33) of Ref. [11] and 
reads 

A « 4mai « i^-. (64) 
2m + + 1 2m + z/ 

Comparing (|63p with (|64|) and noting that n2<Z2 > ??iai, we can see that the presence of the coating 
increases the wavelength of the OSS associated for each mode number. Alternatively, it increases the 
value of the mode number for an OSS of given wavelength. 

We conclude this section by pointing out that we can perform a perturbative calculation of OSS 
in a way that avoids the explicit appearance of the mode number m. According to (I50p and (|51[) . 
tan[4-7rn2(a2 — ai)/A] = tan(7r£). In light of (l4"8j) we can easily compute tan(7r£) and express the latter 
equation in the form 

{( n 2 + Vi + K i) cos(2a) — (rt| — rjx — k\) cosh(2/3)} sin[4-7rrz2A~ (a2 — a±)] 

+2n 2 {r]i sin(2a) + k\ sinh(2/3)} cos[47m 2 A _1 (a2 - a{)] = 0. (65) 



In summary, we have obtained the real equations (j45J) and (J65J) by imposing the complex equation ([34 
We derived the latter by performing first-order perturbation theory on (|25p that determined OSS. It is 
important to note that every solution of (|34p is a solution of (|45p and (|65p . but the converse may not be 
true. We have checked for some concrete examples and found that indeed this is the case. Therefore, 
we solve (|45p and (165|) by fixing all but two of the real parameters entering in these equations and then 
eliminate the solutions that violate (l34l). 



5 OSS of a Concrete Two-Layer Spherical Medium 

In general the refractive index n of an optically active medium depends on the properties of the medium 
and the wavelength of the propagating EM wave. For example, for a gain medium that is obtained by 
doping a host medium of refraction index no and modeled by a two-level atomic system with lower 
and upper level population densities Ni and N u , resonance frequency uq, and damping coefficient 7, it 
satisfies the dispersion relation [8] : 

n 2 =n 2 - „ 2 p (66) 
or — 1 + i^oj 



where Cj := oj/ujo, 7 := j/ujo, oj p := (JVj — N u )e 2 /(m e £oU!Q), e is electron's charge, and m e is its mass 
4 



We can express ui in terms of the imaginary part kq of n at the resonance wavelength Aq := 2ttc/ujq 



according to [8] 

Lb 2 w 2n jK , (67) 

where the approximation symbol means that we neglect quadratic and higher order terms in kq. 
Inserting (I67p in (|66p and using n = rj + in, we obtain 

77 w n + k /i(w), « « Koh(ti), (68) 
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where 



7(1 



Co 2 ) 



j 2 u 



(1 - uj 2 ) 2 + 7 2 w 2 ' Jzy ' ' (l-w 2 ) 2 + 7 2 w 2 ' 



(69) 



We also note that the gain coefficient of such a medium is given by g = —4ttk/X. In particular, we 
can use this relation to express Ko in terms of the gain coefficient go at the resonance wavelength Ao 
according to 

Ao5o 



Kq 



4tt 



(70) 



In the following, we employ (|68p - (|70p to parameterize the refractive indices n i and n 2 that enter in 
the description of our two-layer spherical model. 

In order to explore the effect of the outer shell on the behavior of spectral singularities, we consider 
a spherical dye laser medium confined in a thin spherical shell of higher-refractive index glass. We 
suppose that the refractive index of the glass takes a constant real value and parameterize the location 
of the spectral singularities using the resonance gain coefficient, go, of the dye and the wavelength A. 

Consider confining a Rose Bengal-DMSO (Dimethyl sulfoxide) solution with characteristics [18^ [T9] 



no = 1.479, A = 549nm, 7 = 0.062, g < 5 cm" 



(71) 



in a spherical glass shell of outer radius a>2 = 1 mm, thickness 02 — a\ = 5 /im, and refractive index 
112 = 2.5H Figure [2] and Table [T] show the results of our numerical calculation of the location of 
spectral singularities that use (|25j) . Tabled] also gives the results of our perturbative calculations and 
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Figure 2: Spectral singularities of the Rose Bengal-DMSO dye gain medium (|7ip confined in a spherical 
glass shell of outer radius 1 mm, thickness 5 /im, and refractive index 2.5. The minimum corresponds to 
A = A^ « 549.459 nm that is larger than the resonance wavelength Ao = 549 nm. The grey horizontal 
line represents the experimental upper bound on go. 



demonstrates their good agreement with the numerical results. 

Let <7q^ denote the smallest value of the gain coefficient go that is capable of producing an OSS, 
and A^ 1 ) be the wavelength of this OSS. We recall from Ref. [11] that in the absence of coating A^ 1 ) is 
essentially identical with the resonance wavelength Ao • According to Figure [2] and Table [TJ the presence 
of the glass coating causes A^ to be slightly red-shifted. 

Figure [3] shows a logarithmic plot of the reflection coefficient as a function of the wavelength for 
the sample considered in Figure [2] and Table [T] with go = g^ = 4.85614735 cm -1 . In this case, there 
is an OSS at A = A( 3 ) = 549.56092702 nm that corresponds to the central peak in Figure El Unlike 

3 For the details of developing such a high-refractive index glass, see [20U21| . 
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4.8647 


549.663054 


549.663056 
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9104 


4.8664 


549.152785 


549.152787 


7 


9094 


4.8779 


549.765220 


549.765222 



Table 1: The mode number m, gain coefficient go, and wavelength for the spectral singularities of 
the Rose Bengal-DMSO dye gain medium (|7ip confined in a spherical glass shell of outer radius 1 mm, 
thickness 5 /im, and refractive index 2.5. The subscripts "exact" and "pert." refer to the results of 
numerical and perturbative calculations, respectively. These calculations give the same values for go up 
to six significant figures. 



the other peaks shown in this figure, the hight of the central peak increases indefinitely as we use more 
and more accurate numerical values for the parameters of the system. This is a clear indication that it 
corresponds to a spectral singularity. 



c 
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549.4 549.5 549.6 549.7 549.8 
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Figure 3: (Color online) Graph of reflection coefficient for the outer radius of 1 mm and thickness of 5 
/xm when the gain coefficient is g = 4.85614735 cm" 1 . The central peak represents an OSS. 



Next, we explore the effect of changing the thickness of coating ai — a\ on the A*- 1 -* while keeping the 
inner radius fixed. Figure U] shows the graph of A^ as a function of thickness. As seen from this figure, 
as we increase the thickness, A^ 1 ) undergos an infinite set of jumps that oscillate about the resonance 
wavelength Ao with a decreasing amplitude. Our numerical results show that unlike the wavelength A^ 1 ) , 
the gain coefficient g^ 1 does not experience a noticeable change due to an increase in the thickness. 

Figure shows the location of OSS for a spherical dye gain medium without a coating. Comparing 
this figure with Figure [21 we see that except for the shift in the value of A^ 1 ) , the distribution of the 
OSS in the go-X plane does not seem to get affected by the presence of the coating glass. Exploring a 
wider range of values of A and go reveals a different picture. Figures [6] and [7] show the location of OSS 
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Figure 4: (Color online) Graph of the wavelength A^ as a function of the thickness of coating. The 
inner radius is held at ai = 1.6 mm. As we change the thickness, A^ 1 ) undergoes an infinite set of jumps 
that oscillate about Aq with a decreasing amplitude. 



in a wider spectral range for a coated spherical Rose Bengal-DMSO sample with specifics (171 h . inner 
radius 1.5 mm, and two different values of the coating thickness, namely 5 /jm and 10 /mi. As these 
figures show, OSS are located on curves with multiple local minima. The number of these minima that 
fulfill the experimental upper bound of go < 5 cm" 1 is an increasing function of the thickness of the 
coating. For 02 — a\ = 5 fim (Figure [6]), this number is three, i.e., there are three distinct groups of OSS 
satisfying the bound: go < 5 cm -1 . These appear in the spectral ranges 536.941436-539.544783 nm, 
546.335226-552.571190 nm, and 560.400334-561.815048 nm and respectively contain 41, 93, and 21 
members. Their central member that corresponds to the smallest gain coefficient appear at (A, go) = 
(538.239847,4.601744), (549.435560,3.218290), and (561.106835,4.895167) in (nm,™" 1 ) units. For a 2 - 
ai = 10 /jm (Figure [7J), there are four groups of OSS satisfying the bound on go- They include 34, 47, 46, 
and 30 members with wavelengths ranging over 539.621948-541.782279 nm, 544.992898-548.067753 nm, 
550.878022-553.951299 nm, and 557.516880-559.542755 nm, respectively. The central members of these 
four groups that correspond to the local minima of go have (A, go) values: (540.666848,4.010267), 
(546.525852,3.287748), (552.444850,3.348007) and (558.493430,4.217539) in (mmcm" 1 ) units. In 
particular, the presence of the glass coating not only allows for generating OSS in a smaller gain 
medium, but it produces spectral gaps in the spectral range within which these OSS are located. These 
remarkable observations should in principle be verifiable experimentally. 

Next, we study the effect of changing the radius of the inner core on the location of OSS for a fixed 
value of the coating thickness. Figure [8] shows the graph of the minimum gain coefficient go necessary 
for generating an OSS as a function of a\ for the coated spherical Rose Bengal-DMSO dye gain medium 
with specifications (|71h and thickness 5 /xm. For a-y < a[ min) 0.96 mm no OSS can be created. 
Recalling that for an uncoated sample a^ mm ^ ~ 3.3 mm, this corresponds to a three-fold decrease in 
the size of the gain medium. As we increase ai starting from the critical value 0.96 mm, the minimum 
gain coefficient necessary for generating an OSS, namely g^ decreases. Both of these observations are 
in agreement with our perturbative results. Moreover, it turns out that for fixed values of the coating 
thickness the wavelength A^ 1 ) is not sensitive to the variations of the inner radius. Figure [9] demonstrates 
this behavior. As our perturbative treatment shows the larger the refractive index of coating is the 
smaller the minimum radius a^ mm ^ gets. For example if we use a glass coating with the same thickness 
(5 /im) but a slightly lower index of refraction, say r&2 = 1.93 as in |20j . the value of a^ mm ^ increases by 
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Figure 5: OSS of an uncoated spherical Rose Bengal-DMSO dye gain medium with specifications (|7ip 
and radius 3.3 mm. The minimum gain coefficient necessary for generating an OSS and the corresponding 
wavelength are respectively 4.9815 cm and 549.008 nm. The horizontal grey line represents the 
experimental upper bound on go. There are 66 OSS complying with this bound. 

a factor of 2 (to about 1.8 mm.) 

We have also investigated the formation of OSS in a two-layer spherical model in which the inner 
core is made of a higher-refractive index glass and the outer shell is filled with a dye gain material. 
Using the same glass and gain material and taking the core radius to be 2 mm, we found that the 
minimum shell thickness required for producing an OSS was about 3.28 mm. This coincides with the 
minimum radius supporting an OSS for a sphere filled with the same gain material. Figure [TOl shows the 
location of OSS for the thickness of 3.3 mm as given by f|25j) . The main difference between this figure 
and Figure [5] is that the number of OSS has doubled. To our knowledge this is the only effect of the 
presence of the glass core. We also find that using a glass core with a lower index of refraction does not 
have a sizable effect on the location and number of OSS. 

6 Concluding Remarks 

A gain medium can emit electromagnetic waves provided that we adjust its parameters so that an optical 
spectral singularity is created. The time-reversal of this phenomenon corresponds to the (coherent 
perfect) absorption of electromagnetic waves. This is the theoretical basis of antilasing. The simplest 
example of a gain medium that is localized in space and is capable of realizing an optical spectral 
singularity is a spherical gain medium. It supports spectral singularities in a radial transverse mode 
provided that its radius exceeds a critical value a^ mm \ 

For the typical dye laser material that we consider in Ref. p], « 3.28 mm. In practice 

maintaining a uniform gain coefficient within a spherical sample of this size can be difficult. The main 
motivation for the current study is the idea that coating a spherical gain medium by a high-refractive 
index material can reduce a^ mm \ We have shown by explicit perturbative and numerical calculations 
that this is actually the case. In particular, coating a spherical active dye medium with a glass of 
thickness 5 /im and index of refractive 2.5 reduces a( min ) by about a factor of 3. Furthermore, we have 
found that the presence of coating leads to a much richer structure as far as the location of spectral 
singularities are concerned. In particular, the coating causes a small shift in the wavelength of the 
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Figure 6: OSS of a coated spherical Rose Bengal-DMSO dye gain medium with specifications (|7ip . inner 
radius 1.5 mm, and coating thickness 5 fj,m. The horizontal grey line represents the experimental upper 
bound on go. There are three groups of OSS below this bound. They appear in the spectral ranges 
536.941436-539.544783 nm, 546.335226-552.571190 nm, and 560.400334-561.815048 nm and contain 41, 
93, and 21 members, respectively. 
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Figure 7: OSS of a coated spherical Rose Bengal-DMSO dye gain medium with specifications (|7ip. 
inner radius 1.5 mm, and coating thickness 10 /im. There are four groups of OSS fulfilling the experi- 
mental upper bound on go (the horizontal grey line.) They appear in the spectral ranges 539.621948- 
541.782279 nm, 544.992898-548.067753 nm, 550.878022-553.951299 nm, and 557.516880-559.542755 nm 
and contain 34, 47, 46, and 30 members, respectively. 




spectral singularity that requires the least amount of gain. More importantly, it produces a clustering 
of spectral singularities into groups separated by sizable spectral gaps. This behavior may find an 
application in producing tunable lasers with rather wide spectral gaps. Another interesting possibility 
is to generate laser pulses in different spectral ranges by periodically altering the gain coefficients of the 
active core. For example, as shown in Figures [6] and (8j increasing the gain coefficient so that it passes 
one or more of the local minima of the OSS curves leads to lasing in two or more spectral ranges that 
are separated by gaps of several nm in widthQ 

If we use a glass spherical core and an active dye outer shell, the thickness of the shell required 

4 For an experimental realization of the OSSs we have examined and their possible applications, one may try to use 
the SHG of a diode pumped Nd-YAG (continuous wave) laser at 532 nm wavelength in the plane-wave configuration. To 
achieve reasonably uniform gain, one probably needs to use multiple pumping. The details and a discussion of alternative 
pumping methods are beyond the scope of the present paper and the expertise of the authors. 
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Figure 8: (Color online) Graph of go as a function of the inner radius a\. The coating thickness is kept 
fixed at 5 fj,m. Dashed line corresponds to the experimental bound go < 5cm . For go = 5 cm -1 the 
smallest value of the inner radius that supports an OSS is about 960 \im. For larger values of a±, there 
are OSS with go < 5 cm -1 . 
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Figure 9: Effect of changing the inner radius on the wavelength A^ for the coating thickness of 5 /xm. 
The displayed dots correspond to increments of the radius a± of 100 fim. For continuous variations of 
the radius these dots fluctuate around the line A^ 1 ) = 549.4 nm with a decreasing amplitude as the gain 
coefficient increases. 



to induce an optical spectral singularity is about the a^ rmn ^ for an uncoated spherical gain medium 
consisting of the same dye. The presence of the glass core does not seem to have a significant effect 
except for the doubling of the number of the spectral singularities. We plan to explore the reasons for 
this phenomenon in a more general context. 

The results we have reported apply for a TEM electromagnetic wave propagating in the radial 
direction. As we show in the appendix, we can extend our investigation to study optical spectral 
singularities for electromagnetic waves in the TE and TM modes. It turns out that the results are 
similar to those for the TEM mode. For example the size of the a min remains essentially the same. 

The emergence of an optical spectral singularity is a common feature of any lasing system. This 
includes typical coated or uncoated spherical lasers that can be realized using much smaller, micron size 
dye samples [20^ 122). The reason is that these lasers involve exciting whispering gallery modes [15]. We 
intend to conduct a through study of optical spectral singularities for whispering gallery modes. 
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Figure 10: OSS for a double-layer sphere consisting of a glass core with refractive index ni = 2.5 
and an outer shell filled with a Rose Bengal-DMSO dye solution with specifics (|7ip . The inner and 
outer radii are taken as a\ = 2 mm and a 2 = 5.3 mm. The grey line represents the experimental 
upper bound on go- There are 136 OSS fulfilling this bound. They are located in the spectral range: 
547.969294-550.025813 nm. The OSS with minimum gain has coordinates A (1) = 549.001882 nm and 
= 4.981501 cm" 1 . 
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Appendix 

In this appendix we extend our analysis of optical spectral singularities to the general TE(£, m) 
and TM(f, m) modes of our two-layer spherical system. 

We begin our treatment by writing the solutions of the Maxwell's equations in TE and TM modes 
of a spherical system as follows [17] . 



Lm 



h{r) 



a M (i, m)V x f e (kr)X em + a E (£, m)gi(kr)X lm 



a M {£,m)fe(kr)Xi m - -a E (£,m)X7 x g e (kr)X em 



(72) 
(73) 



where t and m take integer values in ranges [0, oo) and [—£, £] respectively, a E and om are the coefficients 
of TE and TM modes respectively, Z := yVo/eo is the impedance of the vacuum, i(r) is given in ([3]), 
both fi(kr) and gi(kr) have form 



fe(kr),g e (kr) 



A\j v {k\r) + B\n v {k\r) for r < ai, 
M3v{k 2 r) + B 2 n u (k2r) for a\ < r < a 2 , 
Ash^ikr) + B 3 h^\kr) for r > a 2 , 



(74) 



X im (9,<f>) :-- 



LYe tm , L := —ifx V, Yi )ra are spherical harmonics: 
v (a a\ / 2l + l(l-m)!~ pw 



im<f> 



47T U + m)Y 
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and PF 1 are the associated Legendre functions. Note that Xg_ m fulfil the orthogonality relations: 
/ Xp m ,.Xe m d£l = Sei'dmm' and J Xp m ,.(fx X £rn )dVL = 0, where dU := sin 6 d9d<fi, and being proportion 
to the angular momentum operator, L satisfies I?Yg m = £{£ + l)Y^ m . 

Expressing the electric and magnetic fields in terms of their components in the spherical coordinates, 
we have E = rE r (r, 9, 0) + §Eg(r, 9, (j>) + QE^r, 6, <t>) and H = rH r (r, 6, 4>) + 9H e (r, 9, 0) + fiH^r, 9, 0), 
where 
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Imposing the physical matching conditions at r = ai and r = ci2, we find that the tangential components 
of E and H are continuous along the boundaries, i.e., Ef = E™\ Ef = E° ut , Hf = H° ut , and 

For TE modes, we can express these boundary conditions as 



K E 



where for all p, q = 1,2, 
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(75) 



K 



pi 



je(k q a p ) n e (k q a p ) 
je(k q a p ) h e (k q a p ) 



h { p{ka 2 ) hf\ka 2 ) 
h { p{ka 2 ) hf\ka 2 ) 



(76) 



Note that these are respectively identical with (JTSj) and (fT5|) . i.e., = K. pq and = L. Therefore, 
as far as the study of the spectral singularities are concerned, we obtain similar results except that 
v = y/b/2 is now replaced by £. As we see from (|37p . this does not influence the parameters of the 
system. It only changes the mode numbers associated with spectral singularities. 
Similarly, for TM modes, we find the following set of boundary conditions. 
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(77) 
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where, for all p, q = 1, 2, 

(78) 



As we see the only difference between and K. pq is the appearance of the factor 1/ n q in the second 
row of K^. We have shown by explicit calculation that the presence of these extra 1/njj factors does 
not affect the calculation of spectral singularities except for changing the value of the corresponding 
mode numbers. 
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